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ABSTRACT 

We investigate the impact of AGN feedback on the gravitational lensing properties of 
a sample of galaxy clusters with masses in the range 10 14 - 10 15 M@, using state-of- 
the-art simulations. Adopting a ray-tracing algorithm, we compute the cross-section of 
giant arcs from clusters simulated with dark-matter only physics (DM), dark matter 
plus gas with cooling and star formation (CSF), and dark matter plus gas with cooling, 
star formation and AGN feedback (CSFBH). Once AGN feedback is included, baryonic 
physics boosts the strong lensing cross-section by much less than previously estimated 
using clusters simulated with only cooling and star formation. For a cluster with a 
virial mass of 7.4 x 10 14 M@, inclusion of baryonic physics without feedback can boost 
the cross-section by as much as a factor of 3, in agreement with previous studies, 
whereas once AGN feedback is included this maximal figure falls to a factor of 2 at 
most. Typically, clusters simulated with DM and CSFBH physics have similar cross- 
sections for the production of giant arcs. We also investigate how baryonic physics 
affects the weak lensing properties of the simulated clusters by fitting NFW profiles 
to synthetic weak lensing data sets using a Markov Chain Monte Carlo approach, and 
by performing non-parametric mass reconstructions. Without the inclusion of AGN 
feedback, measured concentration parameters can be much larger than those obtained 
with AGN feedback, which are similar to the dark-matter only case. 

Key words: cosmology: observations - cosmology: theory - galaxies: clusters: general 
- gravitational lensing: numerical 



1 INTRODUCTION 

In recent years, strong gravitational lensing by galaxy clus- 
ters has been used to investigate cosmological problems and 
to improve our understanding of cluster physics. Strong lens- 
ing can manifest itself as highly non-linear distortions of 
light from background galaxies (giant arcs) and observa- 
tional searches have revealed many spectacular examples 
around both optically-selected and X-ray-selected clusters 
(e.g. |Zaritsky fc Gonzalez| ( |2003[ ), |Gladders et al.| ( |2003[ )). 
It has been found that the strong lensing properties of clus- 
ters depend on myriad physics - including the cosmological 
parameters, halo substructure and baryonic physics. Strong 



* Email: jmead@ast.cam.ac.uk 
f Email:ljk@ast. cam. ac.uk 



lensing is hence a good technique to investigate both cos- 
mology and cluster physics. 

The frequency of giant arcs around galaxy clusters is 
sensitive to the cosmological model. Using dark matter sim- 
ulations, |Bartelmann et af] ( |1998| ) (henceforth B98) found 
that in ACDM the number of giant arcs is a factor of 5- 
10 lower than in observations. This disparity arises either 
because our cosmological model is incorrect or, more likely, 
that some aspects of the modeling of the source or lens popu- 
lations are incorrect, for example that we are missing crucial 
physics from our cluster simulations. 

A substantial body of work has investigated this latter 
point, using increasingly sophisticated simulations in order 
to reconcile predictions with observations. There have been 
a number of papers investigating the strong-lensing proper- 



ties of dark matter-only clusters (e.g. Ho & White (20051, 
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Hilbert et al. (2007)). The earliest analytic work focused 



on spherically symmetric cluster models, but it has been 
demonstrated that halo triaxiality has a strong impact on 
arc production e.g. by Oguri et al. (20031, who used semi- 



analytic models, and in a complementary analysis by |Dalal] 
et al. (2004) who analysed a simulated cluster sample. 



Other papers have also considered the impact of bary- 
onic physics, in particular models including gas with cooling 
and star formation, and cluster galaxies. Baryons can have a 
significant impact on the central regions of galaxy clusters, 
and since it is precisely this region that is important in arc 
production, one would expect the abundance of giant arcs 
to be sensitive to the baryonic physics of the cl uster. |Puch- 



wein et al. (2005 I, Wambsganss et al. (20071 and Rozo et al 



( 2008 ) found that the addition of baryons to cluster simula- 
tions increases the lensing cross-section by factors of a few. 
Of the cluster galaxies, only the central cD galaxy makes 
a non-negligible contribution to the strong lensing proper- 



ties, increasing cross-section by up to 50% (Meneghetti et 
aL][2003 i 



The baryonic cluster simulations which have been used 
to date in the prediction of lensing cross-sections suffer from 
significant over-cooling, with much higher stellar fractions 



than observed (see the discussion in e.g. Puchwein, Sijacki 
fc Springel| ( [20081 ), |Balogh et aL] ( |2"00Tj ) and~ |Puchwein et 



al. (20101). In the absence of heating, the central regions of 



clusters cool faster than the outer regions, causing gas to 
flow inward in a quasi-hydrostatic cooling flow; this would 
result in a large quantity of cold gas in the central regions 
of clusters, a high star formation rate, and extremely lu- 
minous central galaxies. However, this is at odds with both 
observations of clusters, which do not show such strong cool- 
ing flows, and also of Brightest Cluster Galaxies (BCGs) 
that tend to have evolved, red, stellar populations. There 
is compelling evidence from observational and theoretical 
work that the AGN feedback processes that accompany the 
accretion of gas onto supermassive black holes - residing in 
most galaxies - provide a heating mechanism that prevents 
over-cooling, bringing physical models of clusters more into 
line with obs ervations (e.g. [Valageas fc Silk (19991; Bower 



et al. (2001); McNamara & Nulsen (20071; McCarthy et al 



(2009); |Croton et al.| ( |2006[ )). 

It is well established that strong lensing cross-sections 
can be used to probe the central density profiles of clusters 
- steeper slopes of inner density profiles enhance arc pro- 
duction capability, e.g. | Meneghetti et ah] ( |2007[ ). |Torri et| 



al. (2004) investigated the effects of cluster mergers on the 



strong lensing cross-section, finding that during a merger the 
cross-section can be amplified by up to an order of magni- 
tude. The cross-section for strong lensing can also be boosted 
by as much as 50% for lower mass clusters by line-of-sight 
structure ( Puchwein fc Hilbert||2009 1. An additional factor 
that can affect the cross-section are the properties of the 
background galaxies themselves, such as their redshift dis- 
tribution, size, shape and clustering, e.g. Wambsganss ct al. 



(20041; Li et al. (20051; Gao et al. (20091. 



Departures from ACDM have also been invoked to rec- 
oncile the observed abundance of arcs with the value of 
g 8 (~ 0.8) consistent with cosmic shear surveys and CMB 
observations from WMAP. For example, Bartclman n, Do^| 
|ran fc Wetterich| ( |2006[ ) consider early dark energy (EDE) 
cosmologies, where there is significant dark energy even at 



very early epochs, unlike the case of ACDM or many sim- 
ple quintessence models. Their analytic predictions suggest 
that non-linear structure formation can be substantially en- 
hanced beyond that of ACDM. IFedeli et all 12008b suggest 



that EDE has a significant impact on the statistics of gi- 
ant arcs, bringing predictions and observations into closer 
agreement. More recen tly, |Francis, Lewis fc Under (20091 
and | Grossi fc Springer] (|2009|) carried out the first studies of 
non-linear structure formation in EDE models using N-body 
simulations, finding its impact to be much more difficult to 
detect observationally than implied by analytic models. 

There has also been interest in the effects of bary- 
onic physics on the weak lensing properties of clusters. Fu- 
ture large surveys, such as the Dark Energy Survey (DES; 
https://www.darkenergysurvey.org/) will measure the mat- 
ter power spectrum with unprecedented precision, also al- 
lowing stringent constraints to be placed on the properties of 
dark energy (Albrecht et al. 120061). Studies by, for example, 



|Rudd et al] |2007| > and |White| ( |2004"| ), have indicated that 
the inclusion of baryonic physics can have a significant effect 
on the weak lensing properties of clusters, and hence would 
significantly alter the predictions for the matter power spec- 
trum relative to those simulations which include only dark 
matter. 

It is essential to understand the impact that baryonic 
physics and AGN feedback might have on the gravitational 
lensing properties of galaxy clusters in order to calibrate con- 
straints on their physical properties, and fully exploit their 
potential as cosmological probes. This is important both for 
the interpretation of existing data, and to prepare for up- 
coming large surveys such as DES. In this paper we seek to 
extend the picture of cluster lensing by investigating the in- 
fluence of AGN feedback on the lensing properties of numeri- 
cally simulated clusters. We use recent state-of-the-art simu- 



lations from Puchwein, Sijacki & Springel (20081 that incor- 



porate the growth of black holes (BHs) and corresponding 
feedback processes using the prescription detailed in |Sijack"i] 
et al. (20071. These are the first cosmological simulations of 



clusters to incorporate AGN feedback self-consistently. We 
use these simulations to demonstrate the impact that real- 
istic modelling of AGN physics can have on cluster lensing 
results. 

Throughout this paper a ACDM cosmology with pa- 
rameters tt m = 0.25, = 0.75 and H = 73 kms _1 Mpc _1 
is adopted. The paper is organized as follows. In Section 
2 the cluster simulations that are used are described. The 
numerical methods for the creation of synthetic arcs using 
strong lensing ray-tracing through the simulated clusters, 
and the generation of weak lensing catalogues are outlined 
in Section 3. The impact of baryonic physics on the strong- 
lensing and weak-lensing properties of the clusters are pre- 
sented in Section 4 and in Section 5 respectively. Finally, in 
Section 6 the significance of our findings in the context of 
previous work is discussed. 



2 THE CLUSTER SIMULATIONS 

We consider the gravitational lensing properties of five 
galaxy clusters at z = 0.2 with virial masses in the range 
1.5 x 10 14 — 7.4 x 10 14 Mq, each with three types of physics: 
dark-matter only (DM), dark matter with gas, cooling and 
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Table 1. The virial masses and radii (physical units) of the 
CSFBH clusters used in this paper at z = 0.2. 



Cluster 


M200 


^200 




(1O 14 M ) 


(Mpc) 


CI 


1.5 


0.88 


C2 


2.1 


1.00 


C3 


4.2 


1.25 


C4 


3.1 


1.13 


C5 


7.4 


1.50 



star formation (CSF) and the latter with added AGN feed- 
back (CSFBH). The mass range was chosen to be large 
enough to investigate any strong mass dependence of lens- 
ing properties on baryonic physics. The virial masses and 
radii of these clusters are summarised in Table [T] The virial 
mass and radius are taken to be equivalent to M200 and 
i~200, defined by the mass enclosed inside the radius within 
which the mass density is 200 times that of the critical den- 
sity of the universe at the redshift of the cluster. Note that 
in Table 1, clusters are designated C1-C5 in order of in- 
creasing A/200 with respect to the background density. The 
cluster sample is taken from |Puchwein, Sijacki fe Springel| 
1 2008 ) , in which the simulations are based on the numeri- 



cal prescription described in Sijacki et al. (20071. Puchwein, 
Sijacki fc Springer] ( 2008 1 selected clusters from the dark 



matter-only Millennium simulation ( |Springel et al.|["2005b[ ), 
and performed a number of re-simulations that incorporated 
baryonic physics with and without AGN feedback. DM clus- 
ters that are used in this paper were not already available at 
as high a resolution as the CSF and CSFBH runs, so these 
were re-simulated using the same process as in |Puchwem^ 
Sijacki & Springei (20081 (summarised in Section 2.11. In 
Section |2.2| we recap on the numerical code, and in Section 
2.3|the black hole model is discussed in more depth. 



2.1 Re-simulating Millennium Clusters 

The re-simulations were performed by selecting the La- 
grangian region of the target halo and then populating it 
with a larger number of lower mass particles, adding small- 
scale power up to the new Nyquist frequency. At larger 
distances outside this high-resolution region, increasingly 
more massive particles were used, ensuring that the grav- 
itational tidal field acting on the high resolution region is 
accurately represented. For the simulations with gas, each 
high-resolution dark matter particle is split into a dark mat- 
ter particle and a gas particle, displacing them by half of the 
original mean inter-particle separation. The centre of mass 
of each pair is fixed. 



2.2 The Numerical Code 



The simulation of the cluster sample in |Puchwein, S ijacki 
& Springei! rt2008h follows the same prescription as in | Si— | 
jacki et all (120071) and readers are referred to this paper for 



more details on the simulation and black hole (BH) model. 
The T ree-SPH code GADGET-3 is used (based on |SpringeT 
( |2005[ >), which includes gravity and non-radiative hydrody- 
namics of dark matter and gas components, and also follows 



radiative cooling and heating of an optically thin plasma 
of hydrogen and helium. A sub- resolution multiphase model 
for star formation and the associated supernova feedback 
was adopted (Springei & Hernquist 2003]) . The growth and 
feedback processes of a population of BHs embedded in the 
simulations were self-consistently followed, and this process 
is described in section [2~3l 



2.3 The Black Hole Model 

In the simulations, BHs are represented as collisionless sink 
particles which can accrete gas from their surroundings. Any 
halo above a threshold mass of 5 x 10 10 h" 1 M Q is seeded 
with a BH of mass 1O 5 /i _1 A/0, if it does not contain a 
BH already. The growth of this BH population through 
accretion is then followed. The accretion rate of the BH 
particles follows the spherically symmetric Bondi-Hoyle- 
Lyttleton model JHoyle & Lyttletonl (|1939|; IBondi & Hoyle 



(19441; Bondi (1952)) 



Mb 



4TvaG 2 M£ li p 
(c 2 + « 2 ) 



3/2 



(1) 



where a is a dimensionless parameter (a = 100 in the simu- 
lations we use), p is the density, c s is the sound speed of the 
gas and v is the velocity of the BH relative to the gas. In 
the model, the upper limit of the BH accretion rate is given 
by the Eddington rate. 

There is one additional process through which BHs may 
gain mass - they are allowed to merge with other BHs that 
lie within the smoothing lengths used to estimate the local 
gas density, and have relative velocities lower than the local 
gas sound speed. 

The model for BH feedback accounts for AGN heating 
at both high and low accretion rates. In the model there are 
two modes of AGN feedback: at high accretion rates it is as- 
sumed that the bulk of AGN heating originates in luminous 
quasar activity, whilst at lower accretion rates, correspond- 
ing typically to lower redshifts, AGN heating proceeds via 
radiatively inefficient feedback in a mostly mechanical form. 
The transition between these two regimes is characterized 
by a parameter Xradio = Af B H/M E dd - for Xradio > 10 -2 we 
are in the "quasar" heating regime, and for Xradio ^ 10 -2 we 
are in the "radio mode" feedback regime, which is modelled 
by injecting bubbles into the host cluster. 

In the quasar heating regime a small fraction of the 
bolometric luminosity is coupled thermally and isotropically 



to the surrounding gas particles (Springei, Di Matteo & 
[Hernquist (2005a)), with an amount given by, 



Bfeed = tfLt = EfCrMBHC 



(2) 



In this equation e/ is the efficiency of thermal coupling 
and e r is the radiative efficiency. We adopt e/ = 0.05 and 
e,- = 0.1 as in previous work, to ensure agreement between 
the simulated and observed Mbh — cr* relation ( [Sijacki et al.| 
( |2007[ ), [DTMatteo et ai] ( |2008| ). 

In the radio mode feedback regime, below Xradio, accre- 
tion intermittently produces AGN jets which blow hot bub- 
bles in the surrounding gas. The cluster simulations that we 
use employ a model of radio feedback in which an AGN- 



4 J. M. G. Mead, L. J. King, D. Sijacki, A. Leonard, E. Puchwein, I. G. McCarthy 




Figure 1. Maps of the convergence for cluster C5 using CSF (top row), CSFBH (middle row) and DM (bottom row) for sources at 
z = 3. We show 3 different projections for each physical prescription, reading left to right - xy, yz and xz. The dimensions of the maps 
are 600 kpc in height and 900 kpc in length. The scale bars on the right indicate the ranges of convergence. 



driven bubble is created if a BH has increased in mass by a 
fraction <5bh = SMbh/Mbh- In this case, the energy content 
of the bubble is, 



E h 



l 5 M 



BH , 



(3) 



where e m = 0.2 is the mechanical feedback efficiency. 
The radius of the bubble is modelled by, 



-Rbub — -Rbub.O 



( Ebub/Eb ub,0 

\Picm/picm,o 



1/5 



(4) 



where -Rbub,o> ^Wb.o an d /Oicm.o are normalization con- 
stants for bubble radius, energy content and ambient den- 
sity. 



3 NUMERICAL METHODS FOR LENSING 
SIMULATIONS 

Surface mass density maps are created by projecting the 
matter on to a two-dimensional grid using an adaptive SPH 
projection algorithm. It is also possible to project using a 
triangular-shaped clouds (TSC) scheme. Where the particle 
density is high (i.e. in the strong-lensing region around the 
centre of the cluster) the two algorithms should produce 
almost identical surface densities. We verified that the choice 



of projection algorithm does not affect the lensing cross- 
sections. 

The physical extents of the clusters are much less 
than the distances between the observer-cluster and cluster- 
sources, so the thin-lens approximation is valid. Each surface 
mass density map is determined on a 600 x 600 grid cover- 
ing a square of 1.5 /i -1 Mpc (comoving) on a side, centred on 
the cluster. To create a map of the convergence (k), the sur- 
face density map, E, is divided by the critical surface mass 
density, 

E C = ^^, (5) 



4ttG D ds D d ' 

where D B , D d and D ds are the observer-source, observer- 
lens and lens-source angular diameter distances. In Fig.[l] 
we show the convergence for cluster C5, with all 3 types of 
physics (CSF, DM, CSFBH) in 3 projections, for sources at 
z = 3. 



3.1 Ray-tracing for Strong Lensing Arc 
Production 

The strong lensing properties of a cluster are sensitive to 
the mass distribution in the inner regions. We now describe 
how synthetic maps of distant galaxies strongly lensed by 
simulated clusters are obtained. The ray-tracing for the sim- 
ulation of strong lensing by a cluster, that follows light rays 



The impact of AG N feedback and baryonic cooling on galaxy clusters as gravitational lenses 5 



between the source and observer via deflection at the cluster, 



is performed by an algorithm used in Leonard et al. ( 2009 1 



The mapping between the source and lens (image) planes is 
given by the 'lens equation', 



P = 6 - a (61) 



(6) 



where P are the coordinates in the source plane, are the 
coordinates in the lens plane, and a are the deflection angles. 
The lensing deflection potential ip (9) is related to n through 
a Poisson equation, 



V 2 V - 



2k. 



We can also relate ip (6) to a via, 
Vtp — a . 



(7) 



(8) 



Therefore it follows that k can be related to a in Fourier 
space, 



(9) 



To simulate strong lensing, equation [9] is used to cal- 
culate the deflection angles due to the cluster (on a grid 
of twice the resolution of the input convergence map). To 
achieve greater numerical accuracy, and to mitigate the in- 
fluence of edge effects, the convergence field is zero-padded. 
The output of this ray-tracing program is a map contain- 
ing the lensed arcs. This map is searched for arcs using the 



algorithm presented in Horesh et al. (20051. The algorithm 
makes calls to SExtractor ( jBertin fc Arnou ts 1996), using 
slightly different detection parameters in each call, to de- 
tect objects with an axis ratio above a given value. The 
arc length-to- width ratio is measured as in [Miralda- Escude] 
( 1993 1. 

The giant arc cross-sections are computed using a 
Monte Carlo approach. The cluster field is populated with 
background galaxies placed at randomly assigned positions, 
with exponential profiles and elliptical shaped isophotes 
drawn as Gaussian random deviates, as described in |Leonard] 
et al. (20091. This galaxy population is placed at uniform 
redshift, z g = 3 - although the redshift distribution of the 
source galaxy population is important for determining the 
absolute cross-section, in this work we are interested in the 
relative cross-sections of clusters simulated with different 
types of physics. This population of galaxies is then ray- 
traced through the simulation, using the methods described 
above. In this work, we focus on arcs with a length-to-width 
ratio r ^ 7.5, making no distinction between tangential and 
radial arcs. To calculate the cross-section, 100 populations 
of background galaxies are simulated, each with a density of 
50 arcmin -2 , and the number of arcs produced in each re- 
alisation is recorded. The cross-section is then computed by 
multiplying the fraction of sources that form giant arcs by 
the area of the cluster field over which those galaxies were 
distributed. An error on the cross-section is estimated by 
repeating this procedure 10 times (in all we simulate 1000 
populations behind each cluster), and computing the stan- 
dard deviations of the cross-sections. 



3.2 Weak Lensing Simulations 

Weak lensing is sensitive to the shear field of the cluster 
on larger scales. The weak lensing simulations that are un- 
dertaken to generate synthetic catalogues of weakly lensed 
distant galaxies behind the simulated clusters are now de- 
scribed. These catalogues are used to obtain NFW fits to the 
total mass profile, and non-parametric mass reconstructions 
of the target clusters. 

Whereas the convergence magnifies background objects 
such as galaxies, the shear stretches them. The complex 
shear (a complex quantity with magnitude corresponding 
to the shear amplitude, and phase equal to twice the shear 
position angle) also depends on the second derivatives of the 
lensing potential, 



7 = 7i + ^72 = (^,ii — "0,22) /2 + iip,-_ 



(10) 



where commas denote partial differentiation. The complex 
reduced shear is obtained from the shear and convergence 
through the equation g = 7/(1 — /t). Using the k maps cor- 
responding to clusters simulated with DM, CSF or CSFBH 
physics, the corresponding maps of 7 and g are obtained, by 
relating k and 7 in Fourier space given Eq.[7]and Eq.|10| 

The magnification is given by the inverse Jacobian de- 
terminant of the lens equation, fj.(0) = [det A (6))] _1 . Eval- 
uating this we obtain, 

[{1 -K) -7 2 J 

The ellipticity of a galaxy, describing its shape and ori- 
entation, is denoted by a complex number with a modulus 
related to the axis ratio r of (1 — r)/(l + r), and phase 
being twice the position angle. In the non-critical regime 
(det A > 0) the weakly lensed image of a distant galaxy has 
a complex ellipticity given by, 



(12) 



l + <7*e s 

where e and e s are the lensed and unlensed complex ellip- 
ticities respectively, and g* is the complex conjugate of the 
reduced shear (Bartclm ann fc Schneider|2001 1. 

Galaxy ellipticities are taken to have a Gaussian prob- 
ability density function (Bartelmann & Schncider| |2001 1, 



exp ( 



7R7 2 s [l-exp(-l/a e 2 s )] 



(13) 



Throughout this work we take cr £ = 0.2. 

In the weak-lensing analysis, the simulated background 
galaxies are at z g = 3.0, with random positions on the sky, 
and with ellipticities drawn from the Gaussian probability 
density function given by Eq.[l3] The source galaxy redshift 
was chosen to be sufficiently high so as to ensure that the low 
mass clusters would have cross-sections significantly above 
zero. These galaxies are then lensed by the cluster according 
to Eq. |12[ using the map of g corresponding to the k map 
for a particular cluster, and taking into account any deple- 
tion due to magnification. The slope of the unlensed num- 
ber counts is taken to be s — 0.5, so that the local lensed 
number density til is related to the local unlensed number 
density n via «l = nfJ, B ~. We use a number density of 150 



arcmin , well within the power of future surveys (Lewis 



|fc King||2006| l, although out of reach of present observing 
capabilities. With future observations we ought to be able 
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Figure 2. Spherically averaged density profiles of cluster C4. The left-hand panel shows the cluster simulated with cooling and star 
formation (CSF) and the right-hand panel shows the same cluster simulated with AGN feedback included (CSFBH). The solid green 
curve represents stellar density, the red dashed curve represents the density of dark matter and the dot-dashed blue curve shows the 
density of gas. 



to probe the differences in the weak-lensing properties of 
clusters in great detail. 



4 STRONG LENSING RESULTS 

In this Section our key strong-lensing result, the effect of 
AGN feedback on strong lensing cross-sections, is presented. 
In Fig.[2] we show the spherically averaged density profiles 
of stars, gas and dark matter in cluster C4, simulated both 
with CSF and CSFBH. The stellar density is much higher 
in the CSF run in the central regions of the cluster - a re- 
sult of the 'overcooling problem'. Also note that the dark 
matter density profile has contracted in the CSF simulation 
relative to CSFBH, due to the comparatively large density 



of baryons in the cluster centre (Blumenthal et al, (19861, 
|Gnedin et al] ( [2004] ). 

Fig. [3] shows the image plane magnification for sources 
at z g = 3 in the inner 0.75 /i -1 Mpc of the most massive 
cluster C5 with CSF, CSFBH and DM physics, calculated 
using Eqn.[ll] from the convergence and shear maps. The 
observed properties of strongly lensed arcs depend on the 
distribution of magnification. 



4.1 The impact of baryonic physics on 
cross- section 

The ratios of the cross-sections were calculated for DM, CSF 
and CSFBH clusters and the results are summarised in Fig 
4] Usin g a source plane definition for arc cross-section, |Rozo| 
et all ([20081 and iPuchwein et all ([20051 found that the 



those of DM clusters. We can see that the results of this pa- 
per are broadly in agreement with this result - for instance 
the CSF cross-sections of C5 exceed those of the DM sim- 
ulation by factors of 2-3 (Fig.[4|. However, on inclusion of 
AGN feedback we see that for all but cluster C5, the CSFBH 
and DM cross-sections are comparable. Even for cluster C5, 
the boost in cross-section above the DM cluster gained by 
adding baryons tempered with AGN feedback is < 2. For 
all clusters the reduction in cross-section between CSF and 
CSFBH is statistically significant. For clusters C2, C4 and 
C5 clearly the cross section has been reduced by many times 
the typical error bar. The yz- projection of C3 is anomalous 
insofar as addition of AGN physics caused no notable reduc- 
tion in cross-section. On closer analysis, C3 appears to be a 
merging cluster with large amounts of substructure close to 
the cluster centre - in this sense, this projection of the cluster 
is 'atypical'. The results of CI perhaps warrant closer statis- 
tical analysis; for this purpose, we can use Student's t-test 
with the null hypothesis that the mean cross-section of the 
two samples of clusters (CSF and CSFBH) are equal. It was 
found that in all cases we could reject the null hypothesis 
with greater than 95% confidence. 

The results can be explained by analyzing the density 



profiles of the clusters. Rozo et al. (20081, amongst others, 



inclusion of baryonic cooling boosted cross-sections of the 
most massive clusters by factors of approximately 2-4 above 



have demonstrated that haloes with steeper inner density 
profiles have higher strong lensing cross-sections. Thus, the 
addition of baryons to clusters is expected to boost strong- 
lensing cross-sections, as gas condenses in the centres of 
these haloes, and so the density profiles are steepened signif- 
icantly relative to the DM-only case (e.g. Fig.[2|. In the same 
way, because the CSF haloes have a greater central concen- 
tration of baryons than the CSFBH haloes, the central den- 
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Figure 3. Maps of the magnification for sources at z = 3 in the inner 0.75 h~ 1 Mpc of the most massive cluster C5 with CSF, CSFBH 
and DM physics from left to right respectively. The magnification is shown on a log scale, with the values of the scale indicated on the 
right panel. The black contour corresponds to magnification of 2. 
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Figure 4. Ratios of strong lensing cross-sections for the 3 orthogonal projections of each cluster. In Fig. 4a, 4b and 4c we show the ratio 
of cross-sections CSFBH/CSF, DM/CSF and DM/CSFBH respectively. On the x-axis, the xy, yz and xz projections of CI through C5 
are denoted (cycling through the projections from the lowest to highest M200 with respect to the background density). 
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sity profile is steeper for a CSF cluster than for a CSFBH 
cluster. It is thus unsurprising that the CSFBH haloes prove 
less efficient strong lenses than CSF haloes. The redshift of 
background galaxies will also have an effect on the results. 
For background galaxies with a relatively high redshift, poor 
clusters with CSF gas physics could have a projected sur- 
face mass density greater than the critical density, whereas 
either (or both) of DM or CSFBH clusters could remain sub- 
critical. This would cause a bias. The results of this paper 
will actually be strengthened through use of a lower source 
redshift - CSF clusters, due to their extremely high cen- 
tral density, will almost always exceed critical density; it is 
CSFBH and DM clusters that could first drop to sub-critical 
levels when lower source redshifts are used - this would re- 
sult in an even greater difference between CSF and CSFBH 
cross-sections. To check this result, the cross-sections of C4 
were investigated using a background population at z g = 1. 
In all cases the ratio between CSFBH and CSF cross-sections 
was reduced relative to the z g = 3 case; for example, for the 
xy- projection, the new ratio was calculated as 0.52. As de- 
scribed above, the reduction will be even more pronounced 
for lower mass clusters. 

In order to illustrate the relative effect of AGN feed- 
back for two different clusters, in Fig. [5] the density profile 
of the CSF simulation divided by the corresponding CSFBH 
density profile is plotted for clusters C2 and C4. In C4, the 
density profile is affected comparatively more by the inclu- 
sion of AGN feedback, both in terms of its absolute value 
and its slope. 

As regards to the mass trend, a much larger sample 
will be required to draw firm conclusions. With the excep- 
tion of C5 (which is discussed below), with increasing clus- 
ter mass the ratio CSFBH/CSF decreases, whilst the ratio 
DM/CSFBH is approximately constant. We can physically 
explain these mass trends as follows. The AGN feedback is 
less effective in cluster C5 than C4, insofar as the CSFBH 
cross-section is a larger fraction of the CSF cross-section, 
and exceeds the DM cross-section by a greater amount. A 
possible explanation lies in the fact that AGN feedback 
does not completely overcome central gas cooling in C5. 
Thus, relatively speaking, the AGN feedback has less of 
an impact in C5 than in C4. The stellar fraction of C5, 
even when simulated with CSFBH, is still larger than ob- 
served in reality (Puchwein, Sijacki & Springcl 20081 - if 
the stellar fraction were reduced in C5 to observed lev- 
els this would almost certainly result in a further reduc- 
tion in the CSFBH cross-section, leading both to a reduc- 
tion in the ratio CSFBH/CSF and an increase in the ratio 
DM/CSFBH. On looking at Fig. g]we find that by reduc- 
ing the ratio CSFBH/CSF and increasing DM/CSFBH we 
could possibly fit C5 to both postulated trends - that the 
ratio CSFBH/CSF falls with increasing cluster mass, and 
the ratio DM/CSFBH is approximately constant. On the 
other hand, the AGN feedback is more effective in cluster 
C4 than in C2 - in C4 the CSFBH cross-section is a smaller 
fraction of the CSF cross-section, whilst in both cases the 
DM and CSFBH cross-sections are comparable. This can be 
explained by looking at Fig. [5]- relatively speaking, the AGN 
feedback has had a greater effect on both the magnitude and 
slope of the central density profile in cluster C4 as opposed 
to cluster C2. 

It is well known that for an individual cluster the abso- 
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Figure 5. The density profile of the CSF simulation divided by 
the density profile of the CSFBH simulation for clusters C2 and 
C4. The result is shown for the inner 60 kpc, as it is the den- 
sity profile in this region that primarily affects strong lensing 
cross-sections. The solid line represents C4 and the dashed line 
corresponds to C2. 



lute cross-section of a cluster will vary significantly depend- 
ing on its projection, and we observe this in our clusters. It 
is also worth pointing out that although a mass trend does 
appear to be observed in our results, effects of triaxiality 
and asymmetry effects in clusters could easily mask the true 
mass trend. As mentioned above, a larger sample of clusters 
would be required to investigate this point in detail. 



5 WEAK LENSING RESULTS 

In this Section, the results of the analysis of the synthetic 
weak lensing catalogues are presented, including NFW fits 
to the data and non-parametric mass maps. 



5.1 NFW profile fits 



COS MOMC (|Lewis fc Br idle 2002) is used to fit NFW pro- 
files (iNavarro et al. 19971 to the simulated weak lensing 



data, starting outside the strong lensing regime. The NFW 
profile can be parameterised with a virial radius r2oo, and a 
concentration parameter c. This allows us to define a scale 
radius r s = r2oo/c Inside the virial radius, the mass density 

of the halo equals 200p c , where p c = 3 ^g^ is the critical 
density of the universe at the redshift of the halo. We may 
then write the density profile of the NFW model as, 



p(r) 



<5c Pc 



(r/r s ) (1 + r/r s y 



(14) 



where the characteristic overdensity of the halo, 5c can 
be written, 



5 C = 



200 



3 In (1 + c) - c/ (1 + c) 



(15) 



Using the results of Bartelmann ( 1996 1 we can write 
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Table 2. Concentrations, c, and virial radii, r200; of clusters from an MCMC analysis of simulated weak lcnsing galaxy catalogues. We 
consider all three simulation prescriptions - cooling and star formation (CSF), cooling and star formation with AGN feedback (CSFBH), 
and dark-matter-only (DM). The values shown are for a single projection only (xy). For each value the errors are calculated using the 
67% likelihood contours. 



Cluster 


CSF 


C 

CSFBH 


DM 


CSF 


T200 (Mpc) 

CSFBH 


DM 


CI 
C2 
C3 
C4 
CS 


9.05+J; 6 7 5 3 

= 91 +0.48 

8 91+ 1 ' 61 
°- yl -1.30 

r 91 +0.67 


q 1+1.80 

°-1.52 
9 81+ 2 ' 19 

2 r>6+°- 98 

^■ UD -0.94 

6 9+ 1 ' 01 
°' a -1.00 

4 75+ - 44 

*■ ,o -0.45 


9 c + 1.76 

= 11+0.81 
°- J " L -0.83 
9 Q + 1.16 
z -°-0.97 

5 4+ 87 
J - -0.82 

n cn+0.39 
°- J -0.37 


O.SOt ,;} 1 

1 02+ 07 

i 97 +0.09 
1 ' z '-0.08 

1 20+ 010 
1 - zu -o.io 

1 61+ 004 
J " DJ --0.09 


94+ 19 

U-^_0.15 

1 01+ 009 

1 ^9+0.35 
1,dz -0.24 

1 20+ 011 

, CO+0.06 


og7 +o.i3 

u - 9 '-0.19 

96+ 012 
U - 9D -0.12 

1 90+O.I9 

1 qq + 0.10 

1 '°°-0.09 

1 49+ 013 
i -* M -0.11 



io ,u r 



10 9 =- 



10 B =- 



10 s - 



10 s r 



10 4 




1CTF 



10 9 - 



10 B =- 



10 s - 



10 5 r 



~ i — i i i i i I r 



10 

j-[/r'k P c ] 



100 



10 

/■[/r'kpc ] 



100 



Figure 6. NFW fits to the total density of cluster C4 simulated both with CSF (left) and CSFBH (right). The green solid curve is the 
total density profile of the cluster and the dashed black curve is the NFW fit. A background galaxy density of 150 arcmin - 2 is used, and 
the best-fit NFW profiles are computed using an MCMC approach. As expected, we see that the NFW profile is a particularly poor fit 
to the CSF cluster, while it fits the CSFBH profile somewhat better. 



down the corresponding convergence and shear profiles, 
which are included in Appendix A. 

The probability distribution for the observed galaxy el- 
lipticities, p e , can be obtained from the intrinisc ellipticity 



distribution p e s (Eq. 131 as follows 

2 e 



Pe(e\g) =p e = (e s (e\g)) 



Pe B S ( e lf))- 



eg* - 1| 



(16) 

Using this, we can compute the log-likelihood function from 
the probability density for each lensed galaxy, p ( (ej), 



■X>P«(<*lfl(*0) 



(17) 



which can be evaluated numerically for a trial mass model 
(with corresponding reduced shear field g(0)) given the 
lensed ellipticities. The above likelihood calculation is imple- 



mented as a module in COSMOMC ( [Lewis fc Br idle 2002), 



using it as a generic sampler to explore parameter space with 
the standard Metropolis-Hastings algorithm ( [Metropolis et 
aT] ( |1953| ); [Hastings| < flQT0| >). 

The concentrations and virial radii resulting from this 
analysis are given in Table [2] In clusters CI and C3, the 
effect of AGN feedback is to significantly reduce the concen- 
tration from the CSF case. In the other clusters, the con- 
centration has still been reduced by the inclusion of AGN 
feedback, but not as significantly (insofar as the error bars 
on c for the two types of physics overlap). In contrast, in 
all the clusters, the concentrations of the DM clusters are 
only just below those of the CSFBH clusters. The high con- 
centrations observed for the CSF case are due to the very 
steep inner density profile of CSF clusters - galaxies closer 
to the cluster centre contribute comparatively more to the 
likelihood and thus the MCMC algorithm attempts to fit an 
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Figure 7. Non-parametric mass reconstructions of the cluster C5 performed using the algorithm of Seitz & Schneider (1996). From left 
to right mass-reconstructions of C5 simulated with CSF, CSFBH and DM are shown, using a background galaxy density of 150 arcmin -2 
and a smoothing length of 0.38 arcmin. The size of the box is 0.75 h~ 1 Mpc on a side. Note that it is almost impossible to distinguish 
the CSFBH and DM simulations using the mass-reconstruction alone. 



NFW to this very steep inner density region which leads to 
inflated concentrations. 

As shown in Fig. [6] as is to be expected, the NFW pro- 
file is a very poor fit to CSF clusters, due to the large central 
density spike. The NFW is a better fit to the CSFBH clus- 
ters but is still unable to reproduce the density profile in 
the inner regions of the cluster. At larger radii, the CSFBH 
and DM profiles converge, at a location dependent on the 
properties of the cluster. For the more massive clusters that 
would in practice be studied individually with weak gravita- 
tional lensing, that the c and r2oo values are comparable for 
the CSFBH and DM simulations bodes well for the study of 
such dark matter profiles using weak lensing. 

In general, we expect the concentration of the halo to 
rise as halo mass decreases, and this general trend is ob- 
served, with the exception of cluster C3. As noted earlier, 
in this cluster several dark matter clumps are merging at 
the cluster centre, and this will reduce the concentration 
measured using a weak-lensing analysis when fitting a sin- 
gle NFW profile to the entire mass profile. 



5.2 Non-parametric Mass Reconstructions 

The NFW models that we fit to the shear data above are 
circularly symmetric, and it is clear from Fig. [I] that this 
assumption of symmetry does not hold for many of the pro- 
jections. It is also possible to use the simulated weak lens- 
ing data to perform non-parametric mass reconstructions of 
the clusters, which provides complementary information on 
the detailed structure of the clusters. Non-parametric mass 
reconstructions use a statistical analysis of many weakly 
lensed galaxies to identify over-densities in the mass distri- 
bution. In this paper, we perform finite-field reconstructions 



using the algorithm of Seitz & Schneider ( 1996) 



In Fig. [7] we show mass reconstructions of C5 simulated 
with CSF, CSFBH and DM. From the Figure we see that 
the CSFBH and DM clusters are virtually indistinguishable 
in the mass reconstructions. The CSF cluster shows a higher 
central density peak in the reconstruction, as expected. 



6 DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the impact of baryonic 
physics, in particular AGN feedback, on the strong-lensing 



cross-section of a sample of five simulated clusters. We have 
also briefly considered their weak-lensing properties. 

Using the cluster simulations of |Puchw ein, Sijacki & 
|Spring cl (2008), the first high-resolution cosmological sim- 
ulations to include a self-consistent prescription of BH for- 
mation and AGN feedback, we have simulated gravitational 
lensing by clusters with three different physical prescrip- 
tions, 

• Dark-matter only (DM) 

• Cooling and star formation (CSF) 

• The above plus AGN feedback (CSFBH) 

We have demonstrated that the inclusion of AGN 
physics can have a large impact on the lensing properties 
of clusters by both reducing the strong-lensing cross-section 
relative to CSF physics, and altering the concentrations 
measured in a weak-lensing analysis. The principle conclu- 
sions of this paper are as follows: 

• Inclusion of baryonic physics can have a dramatic im- 
pact on the strong lensing properties of simulated clusters. 
In CSF simulations we find that the boost in cross-section 
relative to DM clusters can be significant - our most mas- 
sive cluster is boosted by a factor of 2-3. The boost is caused 
by the steepening of the density profile, due to an excess of 
stars at the cluster centre. 

• We have demonstrated that with the inclusion of AGN 
feedback, the boost in the cross-section over DM clusters is 
significantly attenuated relative to those simulations which 
include only cooling and star formation (CSF). For all but 
the most massive cluster (C5), where it should be noted 
that feedback is likely under-predicted, the CSFBH and DM 
cross-sections are similar. Even for cluster C5, the CSFBH 
cross-section is < 2 times greater than the DM cross-section. 

• The AGN feedback suppresses the build-up of stars 
and gas in the central regions of clusters, ensuring that the 
CSFBH inner density slope is shallower than in CSF haloes 
- this leads to a reduction in the cross-section. The results 
show some indication for mass dependence of the effect - 
in more massive clusters the CSF cross-section exceeds the 
CSFBH cross-section by more, whilst the ratio of CSFBH 
to DM is approximately constant. To fully investigate this 
interesting result, a larger sample of clusters is required. 

• We have shown that baryonic physics can have a mea- 
surable effect on weak lensing results. The NFW profile is, 
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as expected, a poor fit to CSF clusters and is not an ideal fit 
to CSFBH clusters either. Due to the large central density 
spike of the CSF clusters, measured concentrations can be 
higher than for the CSFBH and DM clusters. In the mass 
reconstructions a similar effect is encountered - we see ev- 
idence for a prominent central density 'spike' in the CSF 
case. 

This work has direct implications for the use of clus- 
ters in cosmological parameter estimation. Arc abundance 
is highly sensitive to baryonic physics, making small sam- 
ples of clusters difficult to use as a tool for precision cos- 
mology. Rather, for small samples, using lensing to probe 
cluster properties is a much more fruitful endeavour. In this 
paper we have demonstrated the sensitivity of cluster lens- 
ing properties to the baryonic physics of the cluster, and so 
if we wish to compare observed arc abundances to those in 
simulations, we must ensure that the baryons are simulated 
as accurately as possible. This has important consequences 
for the interpretation of data from large upcoming surveys. 
The observational prospects for studies of galaxy clusters 
are bright: for example DES will contain ~ 20, 000 clusters 
with mass in excess of 2 x 1O 14 M0, up to ~ 1000 of which 
could contain strongly lensed arcs. This will also allow for 
much larger, homogeneously selected, samples of arcs to be 
obtained. 

Aside from the giant arc cross-section considered here, 
as can be seen from Fig.[3] the DM-only clusters show less 
substructure than the clusters with baryons, in particular 
those without AGN feedback. This will influence the ap- 
pearance of giant arcs that are formed, but also increase the 
cross-section for the production of smaller separation mul- 
tiple images by structure outside the critical region of the 
cluster itself ( |King|2007 1. 

Our key finding is that with the inclusion of AGN feed- 
back, the addition of baryons to dark-matter haloes boosts 
the cross-section for the production of arcs by less than 
previously thought. Clusters simulated with AGN feedback 
have strong lensing cross-sections that are comparable to 
those simulated with dark matter only. The inclusion of 
baryonic physics alone cannot be invoked as a means to in- 
crease the arc abundance much beyond that predicted in 
DM simulations. However, baryonic physics is just one part 
of a larger physical picture, and a realistic comparison with 
observed arc abundances would require not only simulation 
of AGN feedback, but also a realistic source population and 
the inclusion of large scale structure. This is a topic for fu- 
ture work. 
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APPENDIX A: THE CONVERGENCE AND 
SHEAR OF THE NFW PROFILE 

The convergence can be written in terms of a dimensionless 
radial coordinate x — r/r s thus, 



k(x) = K k f(x) 
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The shear is given by, 
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